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1. Introduction 



Simulations in the e regime are complementary to standard large volume simulations. They 
allow to extract low energy constants of the chiral Lagrangian, in some cases with less contamina- 
tions from chiral logs coming from higher order corrections. For a long time it has been thought 
that simulations in the e regime are restricted to chirally invariant lattice formulations. In ref. [jl|] we 
have argued that actually this is not the case, and that simulations in the £ regime can be performed 
also with non chirally invariant lattice actions such as Wilson like fermions. 

In particular in [{[]] we suggested that a combination of algorithmic and theoretical understand- 
ing of Wilson twisted mass makes it possible to actually perform simulations in the e regime with 
Wilson twisted mass fermions. 

Recently it has been shown that with suitable and related algorithmic ideas [Q] it possible to 
reach or get close to the e regime also with improved Wilson fermions. At this lattice conference 
further results in this directions have been presented [Q]. 

In this proceeding we consider a second lattice spacing and we extend to NLO the analysis 
performed in [jlj|. Our setup is a L 3 x T euclidean lattice with spacing a. The lattice action 



S[x,X,U]=S G [U]+S F [x,X,U], 



(1.1) 



is the sum of the so called tree-level improved Symanzik gauge action 
S G [U] = | Y,\ b a L KeTr[l-p( lxl )(jc;/i,v)] +b Y £ Me Tr [l -P^ 2xl \x;ix, v) 



where 



1 

12 : 



with Wilson twisted mass 



Sf[x,X,U] =a 4 ^x(x) Dw + Ws? 3 x(x) 



, (1.2) 



(1.3) 



(1.4) 



where 



1 



£> w = -{7 M (V M + V;)-aV;V M } + m () . (1.5) 

The basic idea of §\§ is that by sampling all topological sectors in the ensemble generation, it 
is not necessary to have an unambiguous definition of topology at finite lattice spacing. To achieve 
this goal it was suggested [|l|] to use a PHMC algorithm ^ treating the low modes exactly and 
reweighting the observables. This could allow to perform simulations at very low quark masses 
without encountering instabilities or metastabilities. 



2. e expansion 

Lowering the quark mass at finite lattice spacing with Wilson-like fermions requires a detailed 
understanding of the interplay between the genuine chiral behaviour induced by the 'pion' dynam- 
ics and the one generated by cutoff effects. A review on the phase diagram and cutoff effects with 
Wilson twisted mass (Wtm) can be found in ref. [^]. In the e regime this is equivalent to saying that 
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it is necessary to understand the coupling of the zero modes with the relevant operators describing 
the effect of the lattice artifacts. The actual values of the lattice spacing, the physical volume and 
the quark mass determine the appropriate power counting, which ought to be used to perform com- 
putations using chiral effective theories. In the continuum the exact integration over the constant 
zero modes can be achieved in the chiral effective theory modifying the p regime power counting, 
in the so called e expansion where the would-be pion mass is small compared to the linear size of 
the box 

-=0(e), Z = 0(e), M K = 0(e 2 ). (2.1) 

As a result of the exact integration the order parameter, or the equivalent ratio R = vanishes 
in the chiral limit at fixed finite volume [g], obtaining restoration of chiral symmetry. One possible 
way to include the effects of the lattice artifacts in this analysis is to include with an appropriate 
power counting the lattice spacing. Here we modify the standard power counting in the following 
way [g] 

M = 0(e 4 ), ~ = 0(e), i = 0(e) a 2 = 0(e 4 ), (2.2) 
where M indicates generically a quark mass. The partition function at leading order reads 

% = J &[U ]e^ Tl i Uo+U o]- C -^ Tr[f/o+f/ t f+ £ f Tr^I/t-f/o)] _ (2J) 

where the scaling variables are 

Zl =c l V = B F 2 m'V, Z2 = c 2 V = - F W ^ , z 3 = c 3 V =B F 2 ii R V. (2.4) 

To argue that this is a proper power counting for actual numerical simulations we list here some 
values 

M~5MeV, a~0.1fm, L~1.5fm (2.5) 
F~90MeV, B ^5.5GeV, |w'| ~ (570MeV) 4 . (2.6) 
Using these values to estimate the relevant scaling variables in this regime one obtains 

MF 2 B V ~ 0.75, a 2 F 2 \w'\V ~ 0.75, 4^ ~ 1, (2.7) 

ar\w | 

which indicates that this is an appropriate power counting. 1 The chiral condensate can be computed 
in the standard way 

r = 1 log ar, zi = o, (2.8) 

and fig. [T] shows the quark mass (left plot) and lattice spacing (right plot) dependence of the chiral 
condensate. We can certainly conclude that the dependence on the quark mass is, as expected, 
smooth, and the cutoff effects are under control. Extension of this computation to NLO including 
standard 2-point functions is currently in progress The power counting introduced is general 
and valid also for plain Wilson fermions (23 = 0). The same power counting could be used to 
develop an expansion with staggered fermions and to check the chiral properties of the spectrum in 
the presence of roots of the staggered determinant. 

1 If the lattice spacing is much smaller a different power counting ought to be used where the lattice artifacts only 
appear at NNLO. 
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Figure 1: Quark mass (left plot) and lattice spacing (right plot) dependence for the single flavour chiral 
condensate normalized with its LO value in the continuum and infinite volume 
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a[fm] 
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6.61(3) 


0.0656(11) 


1.31 


0.00045(12) 



Table 1: Summary of the simulation setup and of the basic ensemble parameters. 



3. Numerical results 

Details of the algorithm used to generate the gauge ensemble can be found in ref. [jl]]. In this 
proceeding we complement the results obtained in [jl]] with a second lattice spacing with a NLO 
analysis. The inversions for the quark propagator have been performed with a stochastic Z2 x Z2 
source located randomly along the euclidean time. Table [I] summarizes the simulation setup. In 
the left plot of fig. g, we show in the first strip the plaquette MC history and its distribution. In 
the second strip we show the MC history and distribution of the lowest eigenvalue, compared with 
the value of the infrared cutoff (horizontal red line) provided by the twisted mass. In the third 
strip we show the MC history of the reweighting factor and its distribution. One crucial parameter 
for stability issues and for controlling discretization errors is the PCAC mass. In the right plot of 
fig. ||, we show the MC history and the distribution of the PCAC mass at xq = 7/4, together with 
the euclidean time dependence of the PCAC mass. It is remarkable that there is almost no sign of 
boundary 0(a) cutoff effects. The analysis gives with the corresponding Z factors [ |l0| ] 

ampcAC =0.00045(12) => aM^(2GeV) = 0.0012(2), (3.1) 

where ^ 

Mf 5 (2GeV) = — M M = ^(Z A m PC Ac) 2 + M q 2 - (3-2) 

We are clearly not at full twist. It is important to remark that this is not so relevant in the regime 
where chiral symmetry is restored. Automatic O(a) improvement [ pT| ] actually holds in a finite 
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Figure 2: Left plot: MC histories and distributions of the plaquette (first strip) smallest eigenvalue (second 
strip) and reweighting factor (third strip). The smallest eigenvalue is compared with the infrared cutoff 
provided by the twisted mass (horizontal line). Right plot: MC history and distribution at xq = T/4 (first 
strip) together with the euclidean time dependence of the PCAC mass (second strip). 



volume and with suitable boundary conditions also for massless Wilson fermions [12]. This is 
somehow related to the fact that in the region where chiral symmetry is restored only O(ampcAc) 
cutoff effects are expected, i.e. very small effects. On the other hand if the mass is of 0(a 2 ) in 
general the latter could be the cutoff effects to become visible. 

3.1 Low energy constants 

The values of the low energy constants (LEC) can be extracted comparing the results of the 
numerical simulations for the euclidean time dependence of basic two-point functions with the 



prediction of #PT [ |13| , |14| ]. In this proceeding we consider the correlation function 

.3 



C p(*o) = 73 LCp(x,x ) 8 ab C P (x,x Q ) = (P a (x,x )P b (0,0)) 
between charged pseudoscalar currents 



(3-3) 



P a {x)=x{x)iy 5 -x{x). 
The euclidean time dependence of the correlation function in #PT is given by 



C P (x ) = a P + -jbp 



y 



where we have defined the following variables 



a? 



-G] (u) 



1 

24 



b P =F 2 B 2 Q 



- ^1 _ I 

y- T 2 ' 



l--Gl(M) 



(3.4) 



(3.5) 



(3.6) 



Details on the definitions of p, u and gi can be found in [13, 14]. We can thus fit the results from 
the numerical simulations with the following fit formulas 



Cp(x ) =A Q +A2y 2 



A 2 

ap = A + — 



2L j 



(3.7) 
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Figure 3: Left plot: Euclidean time dependence of Cp(xo) together with the result of the fit. The solid line 
indicates the fit range, while the dashed line indicates the same curve outside the fit range. Right plot: fit 
results for chiral condensate (first strip) and decay constant (second strip) as a function of the number of 
points included in the fit around the middle point T /2. The circles indicate the actual values quoted in the 
text and the dashed lines indicate the range of values used to determine the systematic error. 

In the left plot of fig. S we show the euclidean time dependence of Cp(xq) together with the result 
of the fit (red curve). The line is solid along the fit range and it becomes dashed where the points 
are not included in the fit. The results of the analysis are 



a 3 L 3 A = (5.94(36)) • 10 



-3 



a 3 L 3 A 2 



(4.81(30)) -10 



-2 



(3.8) 



where the errors have been computed with a nested Jackknife/bootstrap procedure. It is important 
to check the stability of the fit results with respect to the number of data A^data included in the 
fit. This is especially important if we want to make sure that the parabolic time dependence is a 
real feature coming from simulating in the e regime and not just accidental, i.e. coming from the 
standard cosh dependence which can reproduce a fake parabolic behaviour around the T/2. 

In the right plot of fig. ||, we show the stability of the effective chiral condensate and decay 
constant as a function of the number of data points (i.e time slices) around the middle point included 
in the fits. While the chiral condensate shows a remarkably stable result including more points in 
the fit, the decay constant shows a somehow not completely flat dependence on the number of data 
included in the fit. Although this is not worrisome, it might be an indication of a physical volume 
not sufficiently large to suppress higher order corrections. A perfectly well defined way to proceed 
would be to include in the systematic error for F the spread of its value in the region between the 2 
dashed lines. The preliminary result of this analysis is 



roZ 1 / 3 = 0.620(8), r F = 0.220(8)(10) 



(3.9) 



which agrees rather well with results obtained in the e regime using improved Wilson fermions [15]. 



4. Conclusions and outlook 

We are establishing the basic knowledge to simulate with Wilson-like fermions in the e regime. 
To do this we have introduced a power counting to study the e expansion with Wilson-like fermions. 
The LO computation for the chiral condensate confirms the absence of any phase transition, and a 
NLO extension for the condensate and other observables is currently ongoing. 
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Numerical simulations in the e regime with Wtm are feasible using a PHMC with exact 
reweighting. This particular choice allows to lower significantly the quark mass. The extraction of 
LEC such as £ and F becomes then feasible. Moreover the NLO £ expansion is not contaminated 
by chiral logs, which could be a benefit in reducing the systematic errors. 

Computations in this extreme region with Wilson-like fermions require a detailed understand- 
ing of the usual systematics: discretization errors, quark mass and volume dependence. 

We remark that it might be advantageous to combine p and e regime simulations both as a tool 
to attack 2 + 1(+1) simulations, and to further constrain the values of the LEC. 
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